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Abstract 



In order to better understand what to expect from numerical CORE com- 
putations for two-dimensional massless QED (the Schwinger model) we wish 
to obtain some analytic control over the approach to the continuum limit for 
various choices of fermion derivative. To this end we study the Hamiltonian 
formulation of the lattice Schwinger model (i.e., the theory defined on the 
spatial lattice with continuous time) in = gauge. We begin with a dis- 
cussion of the solution of the Hamilton equations of motion in the continuum, 
we then parallel the derivation of the continuum solution within the lattice 
framework for a range of fermion derivatives. The equations of motion for 
the Fourier transform of the lattice charge density operator show explicitly 
why it is a regulated version of this operator which corresponds to the point- 
split operator of the continuum theory and the sense in which the regulated 
lattice operator can be treated as a Bose field. The same formulas explicitly 
exhibit operators whose matrix elements measure the lack of approach to the 
continuum physics. We show that both chirality violating Wilson-type and 
chirality preserving SLAC-type derivatives correctly reproduce the continuum 
theory and show that there is a clear connection between the strong and weak 
coupling limits of a theory based upon a generalized SLAC-type derivative. 
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I. INTRODUCTION 



It was argued in an earlier paper that the Contractor Renormahzation Group(CORE) 
method can be used to map a theory of lattice fermions and gauge fields into an equivalent 
highly frustrated anti-ferromagnet. Although explicit computations were presented only 
for the free fermion theory, it was argued that a corresponding mapping must exist for the 
interacting theory because the space of retained states used for the free theory coincides with 
the set of lowest energy states of the strongly coupled gauge theory. While this argument 
is true, it is obviously important to have a better understanding of the details of how the 
mapping works. In order to get some experience with this process for a theory which is 
well understood we decided to study the lattice Schwinger model (i.e., two-dimensional 
QED), since the exact continuum solution of this model exists. Before diving into the 
CORE computation, however, we first needed to understand the degree to which the lattice 
model exhibits the interesting features of the continuum theory. This paper is devoted to 
an analytical treatment of the lattice Schwinger model with an eye to clarifying the physics 
which underlines the continuum solution and identifying those general features of the model 
which should provide an ultimate check of the correctness of any numerical solution. 

The continuum Schwinger model in addition to being a non-trivial interacting 

theory of fermions and gauge fields, provides a laboratory for studying a wide range of 
interesting phenomena. It exhibits: confinement of the fermionic degrees of freedom and 
the concomitant appearance of a massive boson in the exact spectrum; breaking of chiral 
symmetry through the axial anomaly; screening of external charges and background electric 
fields; infinite degeneracy of the vacuum states of the theory (theta parameters); and the 
ability to produce arbitrary fermionic polarization charge densities by applying an operator 
of the form e^/d^°(^)ig(^) to the vacuum state (due to the anomalous commutator of the 
electric and axial-charge density operators). It is important to ask which of these features can 
be understood in the lattice theory before taking the continuum limit and how complicated 
a CORE computation has to be in order to extract this physics. Although the literature 
contains discussions of various aspects of the model, such as confinement and the axial 
anomaly p|- p!l|] , we are not aware of any systematic discussion of the theory which attempts 
to parallel the derivation of the continuum solution within the lattice framework. This is 
what we do in this paper. 

In order to make the physics as transparent as possible we formulate the Hamiltonian 
version of the theory in = gauge and only then rewrite it within the super-selected 
sector of gauge-invariant states. We then study the Hamilton equations of motion for the 
electric charge density operator, whose form is completely determined by the way in which 
local gauge invariance is introduced into the lattice theory. Obviously, the form of the 
operator equations of motion depends upon the specific lattice fermion derivative and so 
we study this problem for a wide class of different derivatives; in particular, generalizations 



of the so-called SLAC derivative |jT3|, which explicitly maintain the lattice chiral symmetry 



and generalizations of the Wilson derivative |T^, which break the chiral symmetry for non- 
zero momenta. We find that all of these approaches produce a satisfactory treatment of the 
continuum theory, however the detailed physical picture of how things work varies greatly. 

We show that a key issue for connecting the lattice theory to the continuum theory is 
which lattice currents go over to the continuum current operators jo(^) and jo{x). Obviously 
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the local lattice charge density operator, whose form is fixed by the way in which one 
introduces gauge invariance, cannot have this property because the normal ordered version 
of this operator satisfies the identity jo{i)^ = jo{i) for all values of the lattice spacing (since 
only the charges 0, 1, —1 can exist on a single lattice site). On the other hand, as we will 
show, the Fourier transform jo{k) can be treated as a boson operator and the the dynamics 
of the theory tells us that the current operators of the continuum theory are obtained by 
forming an appropriately regulated version of these lattice operators. 

In order to make our discussion essentially self-contained we begin by briefiy reviewing 
the Aq = gauge treatment of the Hamiltonian version of the continuum Schwinger model. 
We discuss: the need for imposing a state condition, such as restricting to gauge-invariant 
states; why only the total Q = sector of the theory can exist at finite energy; and why 
different sectors of gauge-invariant states exist and are labelled by a continuous parameter 
— 1/2 > e < 1/2, which can be identified as a background electric field. Finally, we review 
the Hamiltonian derivation of the fact that the electric charge density is a free massive Bose 
field and the role played by the anomalous commutator of the electric and axial charge 
density operators in the derivation of this result. After reviewing the continuum theory we 
set up and discuss the physics of the lattice version of the Schwinger model in Ao{x) = 
gauge. We then parallel the continuum arguments as closely as possible for a variety of 
fermion derivatives. A careful treatment of the Hamilton equations of motion for the Fourier 
transform of the charge density operator leads to an understanding of how regulated versions 
of these operators go over to the point-split operators of the continuum theory and the sense 
in which these regulated operators can be treated as Bose fields. The difference between the 
way in which things work for generalized SLAC-type derivatives and Wilson-type derivatives 
becomes clear due to this discussion, as does the connection between the strong and weak 
coupling theory for generalized SLAC-type derivatives. 



II. THE CONTINUUM SCHWINGER MODEL 

Hamiltonian formulations of the continuum Schwinger model have been discussed in the 
literature 1^,0] • Our discussion will parallel these discussions to a degree but will differ in 
important details. Our goal is to allow the reader to understand the important features of 
the Schwinger model without unnecessary formalism. 

As we have already noted, the Schwinger model is simply QED in 1 + 1 dimensions, and 
has a Lagrangian density given by: 

C = i^itd,^, + eA,^,)^ - \f,,F>^\ (1) 

In 1 + 1 dimensions there are only three ant i- commuting 7 matrices, 7o,7i,75, and so 
they can be realized in terms of the Pauli a-matrices: 

70 = -iax, 

71 = -iay, 

75 = 7o7i = (^z- (2) 
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In order to enable us to give the most physical treatment of gauge-invariance of the theory 
we choose to work in temporal, or Aq{x) = 0, gauge. Making this choice the Lagrangian 
density becomes: 

£ = iPizd^j^ - eA{x)^,)^P + ^{doAix) f. (3) 

Here, for convenience, we have denoted the spatial component of the vector potential as 
A{x) and dropped its subscript. Eq.(^ tells us that the electric field, 

E{x)=doA{x), (4) 

is the canonical momentum conjugate to A{x) and it has the usual equal-time commutation 
relations with A{x): 

[E{x),A{x')] = ~i6{x-x'). (5) 
Similarly, the fermion operators satisfy the anti-commutation relations 

{iji{x),ijp{x')} = 6{X - x')Sal3- (6) 

It follows immediately that the Hamiltonian in = gauge is 

'E(x? 



H= dx 



+ ip^{x) {idi + ieA[x)) azip{x) 



(7) 



There is an essential piece of the physics of working in = gauge which requires 
discussion. Since we begin by setting Aq = in the Lagrangian, we cannot vary C with 
respect to Aq or OqAq, and so we do not obtain Gauss' law 

G{x) = (d.^E{x) - etl)\x)t^{x)) = (8) 

as an operator equation of motion. In fact, using the canonical commutation relation, Eq.(^), 
we see that 

^-iJdya{y)Aiy) Qf^^^ J dya{y)A{y) ^ (^d^[E{x) + a{x)) - 6^^ 

= G{x) + d^a{x). (9) 
This means that even if we start with a state |0) for which 

G(x)|0) = O, (10) 

we can generate states of the form 

|^^) = e./dM€M(0|^) (11) 

for which 

G{x)\<p^) = dMx)\<Pa)- (12) 
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Fortunately, the operators G{x) (which we identify with the generators of time-independent 
gauge transformations) commute with one another and with H, and so they can all be 
simultaneously diagonalized. Thus we are free to impose Eq.(|l^) as a state condition because 
the Hamiltonian cannot take us out of this sector of the Hilbert space. Actually, we are free 
to impose the more general condition of Eq.(|T2|) for any arbitrary function a{x). What 
all this means is that the canonical quantization of the Schwinger model in = gauge 
produces not one, but rather an infinite number of theories distinguished from one another 
by the fact that they have, in addition to the dynamical fermion fields, different static 
classical background charge distributions p(x)ciass = d^a^x). This shouldn't be a surprise 
because one should be able to formulate QED in the presence of an arbitrary distribution 
of static classical background charges. By quantizing in = gauge all we are doing is 
obtaining all of these possibilities at the same time. 

Since, on physical grounds, we are not interested in formulating the Schwinger model in 
the presence of any non-dynamical charge density, it is customary to limit attention to the 
so-called gauge-invariant states defined by the condition that p(x)ciass = 0. Note that this 
doesn't quite reduce us to a single possibility since all it means is that (9^a(x) = or, in other 
words, a{x) can be an arbitrary constant. If we make such a transformation we shift the 
operators E{x) by a constant, which means that we are free to formulate the theory in the 
presence of a constant background field e. If we worked in finite volume this would amount to 
allowing for the possiblity that there are non- vanishing classical charges on the boundaries; 
i.e. the remaining sectors of the theory differ by a choice of boundary conditions. One key 
question associated with the Schwinger model is whether or not the physics is different for 
different values of the background field. In particular, does the ground-state energy density, 
which is certainly different for the free theory, depend upon the value of e when interacting 
fermions are introduced into the game. 

A simple argument given by Coleman shows that values of e which differ by an integer 
must be equivalent to one another. Before giving the details of the argument it is important 
to note that in one dimension the solution to the equation 

d,E{x) = J2e,6{x-x,) (13) 

j 

for a set of charges ej located at positions xj only has a finite energy solution when the 
total charge J2j^j = 0. This is so because Eq.(|T3|) tells us that in the regions between the 
two charges E{x) is constant and it changes by an amount Cj at each point Xj. If the sum 
of the Cj's is not zero then, assuming the field vanishes to the left of the first charge at xi, 
the field must continue to infinity to the right of the last charge. This means that in order 
to minimize the field energy / E'^{x)/2 one or more of the charges must move off to infinity 
leaving behind a totally neutral system. In particular, if we assume no background field 
then the energy of a pair of particles with charges ±1 separated by a distance s is s/2. In 
the presence of a background field e > the situation is different. When the field is present, 
there is a background energy density equal to If one now separates a pair of charges 

oriented so as to reduce the field between the charges to e — 1, the total change in the energy 
of the system is given by 



where the term —se^ occurs because in the region of length s we have replaced the original 
background field e by e — 1. From Eq.(0) it follows that for e < 1/2 increasing the separation 
between the charges costs energy, while for e > 1/2 separating the charges will gain energy 
(i.e., by moving the charges off to infinity one reduces the background field to e = e — 1 and 
gains an infinite amount of energy). Clearly with that kind of energy gain nothing can stop 
this process from happening and, since the only change in the problem is that now there will 
be pairs of charges at ±cxo, it will continue until the background field is reduced to the region 
—1/2 < e < 1/2. For historical reasons this reduced range of e is usually parametrized by an 
angle 6 = 27re and is one of the two angles which label the exact solutions to the continuum 
Schwinger model 

If we work in the sector of physical states for which 

oo 

Qtot|0) = e J dep(Ol0)=O, (15) 

— oo 

we can solve for E{x) in terms of p{x) 

X oo 

E = el d^piO, Qtot = e J d^piO = 0. (16) 

— oo 

Substituting this into the Hamiltonian we obtain 

H = J dxi!^{x)id^a^i){x) - ^ J dxdyp{x)\x - y\p{y), (17) 

where ■ii){x) = e *-^-°°*^^"^*-^V(2^)- This field transformation enables us to eliminate the term 
A[x) from the Hamiltonian and simultaneously preserve the canonical commutation relations 
of operators iplx), ?/'^(x). It is important to observe that even if we had not been able to 
eliminate E{x) from the Hamiltonian we could have still made this definition but it would 
not have been particularly useful since in that case E{x) would have non-trivial equal time 
commutators with the fermion fields and we couldn't use the canonical quantization rules 
to carry out computations. Note, in what follows we will, by abuse of notation, drop the 
tilde and simply write ipi^) as ip{x). 

The content of the exact solution of this model is that it is the theory of a free boson 
of mass = c^/tt and, moreover, the charge density operator p{x) can be used as an 
interpolating field for this particle because it satisfies a free field equation with the same 
mass. To see how this happens all we need to do is derive the Heisenberg equations of 
motion for p{x). 

The time derivative of p{x) is 

dop{x) = -^[p{x),H]. (18) 

Since p{x) commutes with itself, we use canonical equal time anti-commutation relations for 
the fermionic fields Eq.(^ and obtain: 

dopix) = dj{x), (19) 
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where j{x) = 4'K^)'^z'ip{x). Eg. ([T9|) simply states that divergence of the vector current 
vanishes; i.e., the vector current is conserved. 

The second derivative of the charge density operator is now given by 

d',pix) = -[d,jix),H], (20) 
z 

which evaluates to 

dlp{x) = dlp{x) - ^ J dyidyall/i - 2/2I ip{yi)[-idj (x) , p{y2)] + [-idj{x),p{yi)]p{y2)) ■ 

The key point in the solution of the Schwinger model is the commutator of j{x) and 
p{x'). It is known that this commutator acquires a Schwinger term which we will compute 
by considering Fourier components of the currents: 

p{x) = [ ^e-^^>„ jix) = [ ^e-^'^j,. (21) 
J Ztt J Ztt 

By introducing creation and annihilation operators for the upper Uk and lower dk com- 
ponents of the fermion fields with standard anticommutation relations: 

{ul = 2n6{k - g), {4, d,} = 2Ti6{k - g), (22) 

one obtains: 

f dZ 

[jk, Pi] = J ^ (W-fc '^1+1 - W-fc-g ui-{u-^d)). (23) 

At first sight, this is zero, since the integration momenta / can be shifted / ^ Z — g in the first 
term of the integrand. This, however, is not true. The problem is that the momenta shifts 
can be safely done only in the operators that are normal ordered with respect to the vacuum 
state, otherwise the difference of two infinite c-numbers appears. Since, in this basis, the 
e = Hamiltonian Hq reads: 

Ho = l^k {uluk-dldk), (24) 

the vacuum (the lowest energy eigenstate of Hq) is obtained by filling all negative energy 
states 

|vac) = n4n4|0), (25) 

fc<0 fc>0 

where |0) is the state annihilated by the u^s and d^s. One may see, that for q 7^ —k in 
Eg. (pSD , the right hand side annihilates the vacuum and hence momenta shifts are allowed. 
For k = —q, however, this is not the case, and that can be easily seen by considering 
[jfc, p_fc]|vac). One finally obtains: 

[jk,Pq] = ^2n6{k + q), (26) 
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which translates to: 



Consequently 



[j{x),p{x')] = -dj{x-x'). (27) 

TT 



[-td,j{x),p{x')] = Ul5{x-x'), (28) 



and we obtain: 

dlp{x) = dlp{x) - ^ J dyidy2\yi - y2\p{yi)dl5{y2 - x). (29) 
Integrating by parts twice and using dl\x — x'\ = 2S{x — x'), we obtain: 

d'oP = dip - -p. (30) 

vr 

We see therefore, that the charge density operator p{t, x) satisfies the equation for the free 
field with the mass /z^ = e^/vr. 

Let us take another look at the role of the anomalous commutation relation and the 
gauge invariance in the exact solution of the Schwinger model. First consider the case e = 0. 
The equations of motion 

[pk,Ho] = kjk, [jk,Ho] = kpk, (31) 
allow us to write the free fermion Hamiltonian as a quadratic polynomial in p^ and jk'- 

1 °° 

^0=2!^^ (PfcP-fc + JfcJ-fe) , (32) 



since, combined with the anomalous commutator Eq. , it produces exactly the same 
Heisenberg equations of motion. Since the gauge invariance of the theory allowed us to 
eliminate A{x) from the Hamiltonian once E{x) was replaced by the Coulomb interaction 
written in terms of the operators pk alone, the full Hamiltonian is obtained by adding the 
operator 



00 



^'-Vi^ (33, 



to Hq. Obviously, Hj is also a quadratic polynomial in p^ and therefore, thanks to the 
equations of motion, the anomalous commutator of the spatial and temporal components of 
the vector current and the gauge invariance, the total Hamiltonian is quadratic in pk and 
jk- This makes the theory completely solvable in the continuum. We will discuss just how 
much of this picture survives when one moves from the continuum to lattice version of the 
theory in the next section. 
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To complete the usual bosonization of the theory we observe that pk and jk don't satisfy 
canonical commutation relations, however a simple rescaling remedies this problem and at 
the same time casts the Hamiltonian into a more familiar form. To be precise, since pk has 
no A; = termQ, we can define 



(^k = —j^pk, Ilfc = v^jfc. (34) 
Then, using Eg . (|26|) , we see that 

[nfc,(T,] = 27r5(fc + g), (35) 
and the Hamiltonian takes the form 

H= ^ n,n_, + + -)aka.k . (36) 

Jo ZTI \ 71 J 

Given the canonical commutation relations for 11^ and and this form of the Hamiltonian, 
it is obvious that we are dealing with the theory of a free massive Bose field. 

Let us now turn to the question of the dependence of the theory on the background 
electric field, or rather to the more general question of what happens in the Schwinger model 
if we introduce static classical charges. The remarkable property of the Schwinger model is 
that independent of their magnitude these charges are screened completely. Understanding 
how this occurs will fully answer the question of how the theory depends upon a background 
electric field, since we already noted that having a background field of magnitude — 1/2 < 
e < 1/2 corresponds to having classical charges of magnitude ±e on the boundaries (or 
equivalently at ±oo ). 

From the solution of the theory in terms of j{x) it is easy to understand the screening 
phenomena, since it follows immediately from Eq.(^). Let us consider the Schwinger model 
with two external charges of the opposite sign: 

Pext(a;) = eQext {S{x - xi) - 6{x - X2)) ■ (37) 

As we have seen already, Eq.(^) gets modified to include the external charge density. For 
this reason the part of the Hamiltonian corresponding to the Coulomb interaction acquires 
an additional term and the new equation of motion becomes: 

dip = dip - p^ {p{x) + p^^t{x)) , (38) 

where p'^ = e^/ir. This equation implies that there is now a classical, time-independent 
component of the charge density operator induced by the external charge which satisfies: 

, , 2 ^ f dk cos(k\x — Xi\) , , 
Pind(a;) = -p eQcxt j — _^ ^2 (xi ^ X2). (39) 



^Pq 7^ would imply that the system is not neutral and that would violate the state condition 
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Computing the integral, we obtain for the induced charge density: 

And = -^^^ (e-'^l--il _ e-'^l--^!) , (40) 

which, as advertised, screens the external charge densities. One interesting feature of the 
screening is that two external charges get screened independently from each other . Note 
also that the screening occurs on the scales Ax ~ which for small coupling constant 
can be rather large. Nevertheless, if we now move the external charges off to infinity, so as 
to go over to the sector which in the free theory would have an external background field, we 
see that this field is totally screened in the groundstate of the interacting theory. Moreover, 
since all of the screening takes place within a finite distance of the boundary, there is no 
contribution to the groundstate energy density coming from the background field. 

We should point out that while the previous computation makes it clear that there 
shouldn't be a change in the energy density of the groundstate, it is not at all obvious that 
there isn't a finite change in the energy of the state due to the regions surrounding the 
screened external charge. In fact, there clearly is such a change when the external charges 
are located at a finite distance from one another; however, the question of what happens as 
one moves these charges to plus and minus infinity is a bit subtle. The crux of the issue 
has to do with a definition of the limiting process. As will become apparent in a moment 
the conventional treatment of the Schwinger model amounts to a prescription in which one 
defines the Hamiltonian of the system limit 

H = lim Hn = lim / d^H{^) (41) 

where Q is the closed finite interval Q = [—uj,uj]. With this definition in mind the usual 
prescription is to first take the classical background charges to plus and minus infinity and 
then to take the limit Q oo. Given this prescription it is clear that the total Hamiltonian 
defined in this way never sees the classical screened charges and therefore there is no change 
in the vacuum energy. In order to see that this is the usual prescription which follows from 
bosonization of the model let us go back to Eg. (|5BD and modify it to include the possibility 
of having an arbitrary external classical charge density Pextix). In configuration space we 
obtain: 

H = Jdx Qn(a;)2 + ^{dM^yf + ^Hx) + , (42) 

where e{x) is the function which satisfies the equation 

Pext(a;) = ^<9^e(a;). (43) 



71 



Now, if we set e(x) equal to a constant e we see that all we have to do is define a{x) = a{x) +e 
and the Hamiltonian becomes identical to the one without a background field: 

H = Jdx (\fl{xf + lidxH^))' + ^Hx?) ■ (44) 
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This is the usual way of handhng this issue and so we see that this treatment says that 
the groundstate energy is independent of the external constant background field, which 
corresponds to the prescription we gave above. 

To complete our discussion of the continuum Schwinger model we present another way 
of seeing the screening of the classical background field which doesn't require working with 
the exact solution to the problem, but only the anomalous commutation relation of p{x) 
and j{x). The key to this discussion is the introduction of the conserved gauge-dependent 
current 

j{x)=j{x) + -A{x). (45) 

TT 

Obviously, since A{x) doesn't commute with the gauge- generators G{x) defined in Eq.(|), 
this current mixes states which satisfy different forms of the general state-condition defined 



in Eq.(ll). This means that we should think of j{x) as operating in the full Hilbert space of 
the theory obtained by canonical quantization in Aq = gauge without imposing any gauge 
condition. To show that j{x) is conserved we commute it with the Hamiltonian to obtain 

doj{x) = - \j, h]=- [j{x), H] + - [A{x), H] . (46) 
Now, a slight rewrite of the derivation of Eq.(^) gives 

- [j(x), H] = dojix) = d,p{x) - -d^^p{x) = 9,p(x) - -E{x) (47) 

and since by construction doA{x) = E{x), we obtain 

dojix)-d,pix)=0, (48) 

which means the current is conserved. Integrating this equation over all space we obtain, 
under the usual assumptions about surface terms, that 



dx j{x), H 



0, (49) 



a fact we will use in a moment. 

To understand the significance of the fact that j{x) is conserved imagine that we start 
in a sector of the theory whose lowest energy state satisfies {0\G{x)\0) = 0. Next consider 
the transformed state 

f/(a)|0)=e*/'i«"(«)0-(«)+f^(«))|o). (50) 

We already saw in Eq.(|l^) and Eq.(|l2D that the effect of the term proportional to A{C,) in 
the exponent is to shift the field E{x) so that 

{0\U\a) G{x) U{a)\0) = -d^a{x). (51) 

71 

This equation says that U{a) takes us from a state with no background charge density 
to one with background charge density equal to edx<y{x)/7r. Similarly, it follows from the 
commutations relations of p{x) and j{x) and an integration by parts, that 
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{0\U\a)p{x)Uia)\0) = (52) 

71 

Thus, the total effect of applying U (a) to the vacuum of sector of the theory with no classical 
charges is to map this state into a sector which has a non- vanishing classical charge density 
and at the same time to produce a fermionic charge polarization which cancels it exactly. 
Now imagine that a{x) is chosen as in Fig.|I|. Since dxQ:{x) vanishes except in the two narrow 
regions around xl and xr we see that the effect of this operator is to map the original state 
into one which has equal and opposite classical charge densities around xl and xr and 
induced cancelling fermionic polarization charge densities. As we move xl and xr to minus 
and plus infinity respectively the function a{x) becomes a constant and in the limit, the fact 
that U{a) commutes with H implies that 

{0\U^HU\0) = {0\H\0). (53) 

Hence, the energy of the vacuum of the sector with an arbitrary background field is the same 
as the energy of the vacuum of the sector with no background field, which agrees with the 
previous argument for the bosonized version of the theory. 



III. THE LATTICE SCHWINGER MODEL 

Let us now discuss the Hamiltonian version of the Schwinger model on a lattice. In 
the Hamiltonian formalism time is continuous and space is taken to be an infinite lattice 
whose points are separated by a distance a. As in the continuum, we work in = gauge. 
Furthermore, we introduce fermionic variables ipj^ and ipn associated with each site on the 
spatial lattice and replace the continuous fields A{x) and E{x) by conjugate variables An 
and En associated with the link (n, n + 1) joining the sites n and n + 1. This leads to a 
lattice Hamiltonian of the form 

H = He + Hf, (54) 

where 

/ 

n , v-^n —1 . 

He = -T.EI, H, = Y.ii^irK{n-nUe-^^-'^ ^^(^„0"- (55) 

^ n I 
n,n 

Here the kinetic term K{n — n)af3 is a two- by-two matrix for each value oin — n, the fermion 
fields satisfy the anti-commutation relations 

{{^ir,{^n'f}= 6n,n'Sa,p (56) 

and the link fields satisfy the usual harmonic oscillator commutation relations 

[^n,i?„'] = z5„,„'. (57) 

Note that the fermion fields are dimensionless and in order to make the connection to 
continuum fields we will have to rescale them by a factor of 1/ ^/a to give them dimensions 
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of mass^/^. In direct analogy to the continuum theory, the eigenvalue of the operator En is 
the electric flux carried by the link {n,n + 1). Since, as we have seen, the operator e~*^^" 
shifts the flux on the link (n, n + 1) by e it follows that if we define the normal ordered 
charge density operator to be 

Pn =: i^pliJn) :, (58) 

then the operators 

G{n) = En+i -E„- ep„ (59) 

commute with the Hamiltonian. Hence, similar to the continuum, we are free to impose the 
discrete version of Gauss' law 

G{n) 10) = pf- 10) (60) 

as a general state condition. Therefore we see that the lattice and continuum versions of the 
Schwinger model are essentially the same, in that canonical quantization in Aq = gauge 
gives not one version of two-dimensional QED but rather an infinite number of versions of 
the theory corresponding to quantizing in the presence of an arbitrary classical background 
charge distribution. Note that the form of Gauss' law expressed in Eq. (|60|) requires us to use 
the local charge density operator p„ as the lattice analog of the continuum charge density 
operator. 

Paralleling the discussion of the continuum theory as closely as possible, we focus at- 
tention on the zero charge sector of the space of gauge- invariant states; i.e., the ones that 
satisfy the state condition 

G{n)\(f))=0. (61) 

Once again, in this sector we can explicitly solve for En in terms of p„ and eliminate the 
factors of e*^"^" by incorporating them in the definition of ^jJn■ In this way, in the Q = 
sector of gauge-invariant states, the lattice Hamiltonian can be written as: 

H^H,-'^j:P<r\n-rn\p„. (62) 

n,m 

Because the kinetic term K{n — n)af3 is a function of the difference of n and n we can write 
the Hamiltonian in momentum space as: 

H^Jf^ 4 IZ... + X,.^} i, + ^ 7 I (63) 

—Tv/a —ir/a 



COS ak 



Here we have rewritten the Fourier transform of K{n — n)af3 in terms of two functions 
and Xfc, allowing for a very general class of fermion derivatives. Note that in Eq.(|63D and all 
the equations to follow we have adopted the convention that all momentum space operators 
are normalized in a way that the continuum limit is reproduced by taking a — >■ without 
any additional field renormalization. For example: 
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(64) 



Taking our clue from the discussion of the continuum theory we now turn to the derivation 
of the Heisenberg equations of motion for the current p„. The first step, namely computing 



leads us to identify the result of this computation with the divergence of the spatial compo- 
nent of the vector current (or, alternatively, the time component of the axial-vector current) 
jn- Since the discussion to follow is necessarily a bit detailed it is helpful to summarize what 
it will show us in advance. First, we will see that unlike the charge density operator the 
current jn is intrinsically point-split as a consequence of the equations of motion. Second, as 
in the continuum, the important part of the computation of dojn, by taking its commutator 
with H, involves commuting the jn and p„. This computation will show that one cannot 
solve the lattice Schwinger model exactly for any finite value of the lattice spacing a because 
this lattice commutation relation is not the same as its continuum counterpart. Note that 
this feature is related to the properties of the free lattice Hamiltonian rather than being a 
consequence of the interaction. The same computation will show that the continuum limit 
of the naive commutators does not approach the continuum values for the Schwinger model; 
from this we will see why, on dynamical grounds, one has to study what amounts to a 
point-split version of p„ in order to get the correct physics. 

For the purpose of illustration, let us consider explicit forms of Xk and Zk corresponding 
to a number of popular fermion derivatives. In the case of the naive fermion derivative 
Zk — sm{ka)/a, Xk — 0; in the case of the Wilson fermion derivative Zk — sm{ka)/a, Xk — 
r/a (1 — cos(aA;)); and for the SLAC derivative one has Zk — k, Xk = 0. Given any one 
of these derivatives it is easy to find the one-particle energy levels of the non-interacting 
Hamiltonian Hf by rotating the fields: 



doPn = - 



T [pn, H] , 



(65) 



Xk = Uki^k, 



(66) 



where 




(67) 



^k 

This unitary transformation diagonalizes the Hamiltonian: 





(68) 



— 7r/a 



and if we introduce creation and annihilation operators for the x-fields 




(69) 
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with {ul, Uq} = 2n5{k — q) and {dl, dg} = 2n5{k — g), we obtain: 

ir/a 



Hf= J —Eu {uluu - did,) . (70) 



dk 

—ir/a 

Finally, the vacuum state of the free theory is obtained by filling the negative energy sea; 
i.e., 

|vac)= n 4|0). (71) 

— 7r/a<fc<7r/a 

Given these equations it is a straightforward matter to compute the commutator of H 
with pn to obtain doPn- 

doPn = -[pn,H], (72) 
I 

which in the continuum theory is equal to dxj{x) (where j{x) is identified as the spatial 
component of the vector current, or the time component of the axial- vector current). Com- 
puting the commutator of H with p„ is straightforward but we must say a few words about 
how we identify j„- Basically, in order to maintain the parallel to the continuum discussion 
we define the quantity equal to doPn as the lattice derivative of i.e., 

dopn = - (jn+l - jn) ■ (73) 

With this identification, the algebra of 7 matrices in two dimensions ensures that the spatial 
component of the vector current coincides with the temporal component of the axial current, 
and therefore all the currents we are going to work with appear to be defined. Clearly, 
different lattice fermion derivatives will produce different definitions of the spatial component 
of the vector current operator, an inescapable consequence of the Heisenberg equations of 
motion. 

To derive an explicit form for jn, we Fourier transform Eq. ([73|) . Defining 

n 

we obtain: 

doPk = ^[pk,H]. (75) 
Writing the right hand side of this equation as: 

- E (jn+l - Jn) e*''"" = Jfc, (76) 

a „ a 

defines the Fourier transform of the spatial component of the vector current. Explicit com- 
putation of pk yields: 
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■K ja 



Pk 



27r 271 



(77) 



-ir/a 



where is the lattice 5-function which imphes the momentum conservation modulo 

271 /a. Focusing, for the sake of definiteness, on momenta k > 0, one finds: 



Pk 



dfci 

~2^ 



'^t'^ki+k + 



n/a 



~2^ 



i'l^i'ki+k-2n/a- 



(78) 



-Tr/a 



w/a—k 



It is now a straightforward matter to compute the spatial component of the vector current 
using Eq. ([75|) : 



ae 



ika/2 



3k 



where 



2sin(aA;/2) 



■n/a—k 



— 7r/a 



—iljl^M{ki,k)ipk,+k-+ J —^i^M{ki,k)^k^+k_2^/a. 

■K ja—k 



(79) 



(80) 



and we have used the fact that and are periodic functions with the period 27i/a. 

From the continuum solution of the Schwinger model it is clear that we should focus 
on the Schwinger term appearing in the commutator [j^,Pg], since it is the source of the 
anomalous Heisenberg equation of motion and the reason for the mass of the photon being 
non-zero. As we saw in the previous section it suffices to take the vacuum expectation value 
(vac| [j^, Pq] |vac) in order to compute the Schwinger term. Direct computation yields the 
following result: 



(vac|[jl,P,]|vac) = 27r5(A;-g) W{k), 



(81) 



where the function W is: 



w 



ae 



-ika/2 



2sin(aA;/2) 



(2^fci — ^fei-fc — ^fci+fc) cos6'fc^ + {2Xk^ — Xk^_k — Xfc^+fc) sin6'fc^ 

(82) 

To compare the result of this computation with the continuum result we take the limit 
a — > 0, in which case Eq.(^) simplifies and one obtains: 



k 

hm W = — 

ak^O 71 



d^ 



This equation gives the a — limit of the anomalous commutator for a general lattice 
fermion Hamiltonian and is therefore useful for the analysis of the continuum limit of the 
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various choices for the fermion derivative. To get a feehng for how things work let us consider 
several specific examples. 

Let us begin with the case of the naive lattice fermion derivative, where = sin^, 

= 0, E"^ = I sin^|. In this case we obtain: 



k f 2k 
\imW = /d^sin^ = . (83) 

ak^O IT J n 



This shows that the anomalous commutator is two times larger than the continuum result, 
which implies that in the a ^ limit the mass of the photon is two times larger than in 
the continuum theory. In principle, this result should have been expected since the lattice 
theory with the naive fermion derivative has an exact SU (2) symmetry for all values of a 
and as a consequence of this symmetry the fermion spectrum is doubled as is evident from 
the form of Ek- Thus, it follows that the continuum limit of the naive theory is not the 
original Schwinger model, but rather an 5'?7(2)-Schwinger model which is known to have a 
photon mass which is 2e^/7r. 



In case of Wilson fermions we have = sin^, = r(l — cos^) and E^^ = ^ Z"^ + X|. 
By explicit calculation one finds: 

2k 

\imW = , (84) 

afc^O vr 

Once again the result is two times larger than the continuum oneQ, but in this case the low 
energy spectrum is clearly undoubled and the reason for the discrepancy between the lattice 
and continuum results must be different. 

In order to clarify the underlying physics, it is instructive to consider somewhat un- 
conventional fermion derivatives. Let us begin by considering a modified SLAG fermion 



derivative [|TT|. Consider the free fermion Hamiltonian defined by {k > 0, Z^k = —Zk): 



Zk = ke f/i- -k\ + (^-k]e(k- , = o, (85) 

\ a J jjL — 1 \a J \ a J 

and Ek is equal to \Zk\. A plot of E^ is shown in Fig. 2. Computing the second derivative of 
Zk one obtains: 



^The fact that the anomalous commutator for Wilson fermions is r-independent in the a — > 
hmit is a bit of a miracle and we do not quite understand the reason for that. Note however, that 
a similar situation has been observed in earlier calculations of the chiral anomaly in the lattice 
Schwinger model with Wilson fermions Q. It is generally accepted that the continuum limit of 
the Schwinger model with Wilson fermions gives correct anomaly and correctly reproduces all 
other continuum results. However, there is no direct contradiction with our statement since, as we 
explained and as our result seem to illustrate, different currents lead to different results. 
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and the anomalous commutator becomes: 



limW = — + . (87) 

ofe-^O TT \ 1 - IJ, J 

To understand the information encoded in this form of the anomalous commutator let us 
consider what happens in the continuum limit. It is clear from the plot of for this modified 

SLAC derivative that two species of fermions survive in the limit a — >■ 0. Note however that 
dEk/dk is quite different for the two linear regions of the spectrum, which means that the two 
species propagate with very different speeds. The anomalous commutator is really the sum 
of two contributions: one coming from < k < /nr/a and the other from /nr/a < k < tt/o 
and these contributions can be easily identified with the different fermions. Since both 
fermions are charged, they both contribute to the anomalous commutator and to the mass 
gap. 

Given the simple nature of this fermion derivative it is clear how to separate the contri- 
butions of the two fermion species to the total current. The easiest way to do this is to put 
a sharp momentum cut-off somewhere below and above the turning point //7r/a. With this 
prescription we write the charge density operator as a sum of three contributions: 

Pk=f^+fi^+pf\ (88) 

where 

Pk^ = / 2^^^i^'=i+'=' P^^^ ^ J + / ^V'fciV'fci+fe-27r/a, (89) 

—\x\'Kla |fci|>/i27r/a i^ja—k 

and //i < // < //2- The p\ provides for the remaining contribution to the charge density 
and it is only sensitive to fermions with momenta /xivr/a < < /X27r/a. Now, following our 
previous argument, we define the corresponding spatial components of the vector current by 
explicitly commuting the above charge densities with the Hamiltonian. 

A straightforward computation shows that in the limit of vanishingly small lattice spacing 
the anomalous commutators of the above currents are given by: 



^•«V,pW] = -^27r<5(^-,), 

TT 



[(jf)''/^f^] = o, 

\(2f)\pf\--c\2^K^-^\ (90) 

where we introduced c = — /i), which is the velocity of the fermions in the region 
A; ~ TT. Note that in all of these formulas there are also non- vanishing normal-ordered 
operators coming from large momentum excitations which we have not displayed. These 
operators annihilate the vacuum and one can argue that they are unimportant for small k 
physics. This final point, which is intimately related to the sense in which jk can be treated 
as a boson operator, merits elaboration and we will return to it immediately after completing 
our discussion of the equations of motion. 
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Proceeding with our computation of the equations of motion for p^^"* and p^^"* we obtain: 



vr 

-9oVf = c^A:Vmc^pr, (91) 

TT 



where the total charge density operator appears on the right hand side of these equations and 

(2) 

SO pI is still included. Consistent with the point made above and subject to the discussion 
to follow we will set it to zero, since for small energies fermions with such high momenta are 
not excited. 

(3) 

From the equations for the commutators we see that the field p^ is not canonically 
normalized and so we introduce the new field p\ = 1/ \^p^j^^ and its commutation relation 
with the corresponding current becomes canonical. The equations of motion become: 



5oVt^ = c^A:^Pr + -(v/^pincpf)- (92) 



The energies of elementary excitations can be determined from the eigenvalues of the matrix 

j^_(k^ + eVvr e^TT \ . 

■^^~[ eVvr c^k' + c e^ji: ) ' ^^^> 

The matrix is easily analyzed in the limit of large c which corresponds to the large slope of 
the fermion derivative in the region ~ vr. Note that the original SLAC fermion derivative 
corresponds to c — oo limit. 

It is then easy to see that there are two different limits in this equation. For small 
momenta ck^ <^ e^/vr, there are two eigenvalues: 



El = ^ck, E2 = VTT 



vr 



Hence there is a zero mass eigenstate which is the Goldstone boson of the theory. The other 
excitation is the massive one. If we consider the limit c — 00, the region of momenta sensitive 
to the Goldstone mode shrinks to zero (see Fig. 3) and the mass of the other excitation goes 
to infinity. 

For larger momenta, ck^ ^ e^/vr, the mixing of two states becomes small, they propagate 
independently and their energies are given by: 




e2 



^i = Y^' + 7' E2 = ^c^k^ + c-. (94) 

In this momentum region the energy of the lower excitation approaches the result of the con- 
tinuum theory of a bosonic field with the mass e^/vr, the other excitation becomes infinitely 
heavy and decouples explicitly (see Ref. 0]). Hence, as we approach the limit c ^ 00 
(which is equivalent to the orginal form of the SLAC derivative) , the continuum limit of the 
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theory has a massive boson of mass = e^/vr and an isolated state at = which can be 
identified with a seized Goldstone mode [Q. 

Since the main purpose of this paper is to provide an analytic framework for CORE 
computations to follow we should point out that the fact that there is a Goldstone mode 
when c^/tt ^ cA;^ is quite significant since it relates to the whole question of how the strong- 
coupling limit of the model connects to the weak coupling theory and whether one can get 
the correct physics by projecting onto the sector spanned by the e — > cxd eigenstates. We 
will have more to say about this point in the conclusions, but first we should complete our 
discussion of terms we ignored in the commutator of \jkipq\- 

While this preceding argument leading to Eq. (p^) makes the lattice discussion look re- 
markably like the continuum Schwinger model, we are not really finished. The issue which 
still needs discussion relates to the interpretation of pk and jk as boson fields. This is more 
than an academic issue. Although, as we have shown, the vacuum expectation value of 
the commutator of jk and Pq gives the required Schwinger term, computation of the full 
commutator contains an extra piece which, if it has non-vanishing matrix elements between 
states whose energy remains finite in the a — limit, ruins the interpretation of jk and Pq 
as boson fields. This is the issue which we will now address. 

The commutation relation for j^. and Pg, which is valid for arbitrary lattice spacing, reads: 



where the normal ordered operator 0(/c, q) has the form: 



> + {k q). 



(96) 



We will now argue that even though the term 0{k,q) is not explicitly suppressed by 
a power of a, nevertheless this operator does not contribute to the dynamics of any state 
whose energy remains finite as a — 0; in particular, any state which can be created by 
applying arbitrary powers of pk to the groundstate of the theory. As we pointed out in the 
Introduction, the operator p„ cannot be considered a boson operator since p^ = pn, thus 
arbitrary powers of pn can produce at most three linearly independent states when they are 
applied to the groundstate. The situation is quite different for p^^^ and p^^^ which are sums 
of pnS and therefore, for an infinite lattice, will not satisfy an identity of this type. 

The argument begins by considering the non-interacting theory and thinking of 0{k,q) 
as acting on a Hilbert space constructed by applying polynomials in the current operators to 
the groundstate of the free theory. For values of k and q which are small compared to vr/a, 
the operator 0{k, q) can only act on the part of a state which contains left and right moving 
fermions with momenta k ±fii7T/a since it has to first absorb a high momentum fermion 
and create another one with a momentum which differs by a small amount. This means that 
in order to have a matrix element of this operator between states generated by polynomials 
in and p^^"* these states have to have non-vanishing components having fermions with 
high momenta. Thus, we have to ask how such components can be generated? 
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Both p\. and p). are bilinears in fermion creation and annihilation operators and, being 
normal ordered, can only absorb a fermion at one momentum and create a replacement at 
another momentum. Generically these operators are of the general form: 

= / ^ + 4i4i+fc) (97) 

and since for the modified SLAG derivative = 0, the vacuum, as in the continuum, is 
given by Eq.(^). If we now, for the sake of definiteness, consider 

pSo|vac)=p« n4n4|0) (98) 

fc<0 fc>0 

(i) 

we see that almost all terms in annihilate the vacuum state. The only terms which act 
non-trivially are ones were either u^i+fc or d^^^k can absorb a particle and then either or 
4^ can create a particle. Clearly, for small A; > only the d^^ terms can act, since if A;i < 
and ki + k > then dk^+k can absorb a d from the vacuum state and dl_^ can create a d. The 
Uk terms cannot act non-trivially because in order for Mfci+fc to absorb a particle, ki + k has 
to be less than zero, in which case ki < and therefore ul_^ annihilates the resulting state. If, 
however, k < then it is the term which acts non-trivially and the corresponding 

d term annihilates the state. In either event the important point is that the p^*^ only creates 
and absorbs particles from the vacuum which are within a distance \k\ of the top of the 
negative energy sea (i.e., the fermi-surface) thus creating a particle anti-particle pair. 

The next step is to see what happens if we apply p^^ to the state we just generated. 
What we get is 

Pfe^Vac) = J dkidk2 (ul^Uk^+k + dl^dk^+k) (uluk,+k + dl^dk,+k) |vac). (99) 

(i) 

It should be clear that for almost all ki and k2 in the allowed region p^ creates two low 
momentum particle anti-particle pairs and in fact for given allowed ki and k2 there are 2! 
ways of getting the same two-pair state; however, for a given ki there is exactly one value of 
for which one can create a higher energy one-pair state by absorbing one of the particles in 
the pair created by the first application of p^*^ and promoting it to higher momentum. From 
this it follows that the factor needed to normalize this state is greater than 1 / v^. Similarly, 
if one hits this state with another power of p^.*"* almost all of the terms would create three 
low energy particle anti-particle pairs and each of these three pair states would be created 
3! times. As in the previous case however there would be a single term which could promote 
the previous single higher energy one-pair state to yet higher energy. Note, however, that 
since the normalization of this state would have to be larger than l/\/3l (which we are 
beginning to recognize as the normalization factor which goes with a three boson state) 
the coefficient of this higher energy single-pair state appearing in the normalized version 
of the state created by p^^\ is getting smaller each time. If one now imagines carrying 
out this process p-times the argument generalizes in the obvious way. The state obtained 
by applying p powers of p^*"* to |vac) is going to be mostly made of p different low-energy 
particle anti-particle states, each of which will be arrived at in p\ ways. Furthermore, there 
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will be a single particle anti-particle pair state with individual momenta p times larger than 
k. Since now the normalization factor of this state is bigger than p\, the coefficient of this 
single higher energy pair state is getting very small relative to the part of the wavefunction 
made of p low energy pair states. A more careful discussion of this point would also take 
into account the fact that the same procedure will generate two-pair, three-pair, etc., parts 
of the wavefunction. However, the point is that if we keep k < A, where A is a maximal 
energy we wish to consider and Aa ^ in the contimium limit, then in order to achieve the 
fermionic level with momenta Hiir/a, one should create a state jj^|vac) with 

/ilTT /iiTT 

A ~ — — ~ — > oo. (100) 

ak aA 

The energy of this bosonic state is ~ £^typ ~ HiTi/a — > oo; and it is easy to see that the 
probability of finding a single high momentum pair state equals to 1/A^! — > 0. The factors of 
p\ which appear in the normalization of the j'^jvac) states thus produce the explanation of 
both why we can think of p[*^ as a boson operator and why 0{k, q) has no significant matrix 
elements between normalized states generated by applying arbitrary powers of p^*^ to |vac). 

The main point of the above discussion is that as we approach the continuum limit the 
essential physics of the model is taking place near the top of the negative energy sea and so 
it is useful to limit our attention to modified operators p^*^ that only have support in these 
regions. For the model based upon a modified SLAC derivative we saw that, since the low 
energy spectrum of the theory was explicitly doubled, the non-split fermion current really 
was made up of two parts: the first, coming from states near A; ~ and the other from A- ~ tt. 
Leaving aside the complications related to the existence of the Goldstone mode, we see that 
the dynamics of the theory tells us that the current constructed out of the fermionic fields 
with small momenta is essentially the current with the correct continuum limit. The large 
momentum part of the current decouples from the continuum limit after the limit c ^ oo is 
taken. From this point of view, we see that the dynamics of the lattice model tells us that in 
order to take the continuum limit of the theory we have to restrict attention to only a part 
of the unregulated lattice current. This is essentially equivalent to adopting a point-splitting 
procedure for defining the current in the continuum theory. 

Though these peculiarities have been made obvious because of the explicit doubling, 
our calculation of the anomalous commutation relation shows that for the non-point-split 
currents the large momentum modes do not decouple automatically, even without fermion 
doubling. To have explicit decoupling one has to construct the currents by explicitly cutting 
off the region of large momentum. If in the small momentum region the fermion derivative is 
sufficiently continuum-like (i.e., linear) then we are assured that: the current constructed in 
this way will have correct anomalous commutation relations modulo corrections suppressed 
by inverse cut-off; the equations of motion for this current will be free equations of motion 
with an additional source term given by high momentum fermionic modes; if the cut-off is 
sufficiently large in physical energy units (as opposed to lattice units), such a current operator 
will correspond to a continuum bosonic degree of freedom for all low-energy purposes. 

To conclude this section let us consider a simple example of what we will refer to as a 
perfect Wilson model for the fermion derivative. It is defined by 



k9 



TT 

2^ 



k) + ^sin(. - ka) 0{k-£), X,^ ^cos {n - ka) [k - £) , 
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where we have defined and for A; > 0, and assumed that Z^^ — Z^ and X^^ — ^k- 
The one-particle energy spectrum is shown in Fig.4. 

One may easily check (using our general result for the anomalous commutator) that for 
this model, in the limit a — > 0, the commutator for the non-split currents [jk, Pq] differs from 
the continuum limit, even though in this case there is no doubling and the theory remains 
continuum-like up to momenta A; ~ tt/ (2a). Note, however, that if we construct a low energy 
current by restricting to the linear region of the derivative function, as in the case of the 
modified SLAG derivative, we can guarantee that the high momenta modes do not have an 
infiuence on the dynamics of the low energy current and we can verify that this low-energy 
current and its time derivative satisfy the desired anomalous commutation relations. Once 
we establish this fact we can proceed to derive the Hciscnberg equations of motion for this 
low-energy (or regulated) version of the current and make the connection to the continuum 
theory. Conceptually, our example of the perfect Wilson fermion derivative is very close to 
Wilson's original proposal. All we have done is to enlarge the region of momentum space 
in which the lattice derivative looks identical to the continuum derivative so as to make it 
easier to sec why this type of fermion derivative works once the proper current operators 
have been identified. 

IV. CONCLUSIONS 

As we have said in the Introduction our aim is to use the lattice Schwinger model to 
test the idea that one can use CORE methods to map gauge theories into highly frustrated 
spin antiferromagnerst and then use the same methods to study these spin systems. The 
Schwinger model is a very good place to test this notion since the continuum model exhibits 
a rich spectrum of physical phenomena, anomalous commutators, background electric fields, 
charge screening, etc., and so it is important that any numerical treatment of this model 
be able to see these effects. We had several major goals in this paper: first, to get some 
analytic control over the physics of the lattice Schwinger model in order to understand 
which features of the continuum theory we might expect to emerge easily from a numerical 
computation and which might be difficult to obtain; second, to gain a feeling for how much 
of this physics we might hope to see if we first use CORE to map the lattice Schwinger 
model into a highly-frustrated generalized antiferromagnet and then to analyze the physics 
of that spin system before carrying out detailed numerical computations; third, to get a 
better understanding of how the low-energy physics of the lattice system depends upon the 
choice of fermion derivative and why, on dynamical grounds, the lattice currents of interest 
are those which correspond to continuum point-split currents. 

To accomplish our goals we studied Hamiltonian formulations of both the continuum 
and lattice Schwinger model and then, by paralleling the solution of the continuum version 
of the theory in the lattice framework, identified those features of the lattice theory which 
differ from the continuum theory and identified the operators of the lattice theory which go 
over smoothly to their continuum counterparts. It became apparent from the treatment of 
the equations of motion for the various forms of the charge density in the lattice theory that 
getting the right behavior involves showing that one is close enough to the continuum limit 
so that the appropriately defined currents act as bosons; in other words, it is not sufficient 



23 



to only show that there is a gap between the vacuum state and the first excited state and 
that it numerically appears to be of the order of e/y^- At a minimum one should be able 
to show that the operators 0{k,q) have negligible matrix elements between the computed 
low-lying states of the theory. 

A surprising outcome of this work was the fact that almost any fermion derivative works 
for the study of the Schwinger model. As we have seen, the c — > oo limit of the chirality 
conserving modified SLAC derivative and the perfect Wilson derivative had essentially the 
same low-energy behavior. The interesting fact was that the dynamics of the system, while 
different for the two cases, managed to automatically eliminate spurious degrees of freedom. 
Basically this says that we can use any short-range derivative, either chirality preserving 
or chirality violating, to carry out numerical studies of the lattice Schwinger model and by 
comparing them get additional control over how well the numerical methods can be expected 
to converge. 

Finally, and most pertinent to our eventual goal, is the fact that the discussion of the 
modified SLAC derivative shows that the trick of using CORE to map the system into a 
frustrated generalized antiferromagnet will preserve the relevant low energy physics. The 
reason for this is that the CORE method is based on defining the set of retained states to 
be those states which have zero energy in the limit e — >■ oo. This, of course, requires that for 
these states the Coulomb term vanishes. In other words, these are the states for which the 
normal ordered charge density operator p„ is zero identically. (This set of states is generated 
by selecting from the four possible states per site, only the two states having zero charge and 
then taking tensor products of all of these states.) Note that for large these states are all 
degenerate to order and this degeneracy is lifted by the kinetic term which acts on them 
by creating a pair of separated charges and then acting a second time to bring them back 
to a neutral state. A second order degenerate perturbation theory calculation shows that 
the low energy theory in the large limit is that of a Heisenberg anti-ferromagnet, which 
means that in this limit the theory is that of a massless particle. Going back to Eq.(P3D 
we see that perturbing in the kinetic term is the same as taking e^/vr to be much greater 
than k"^ and ck'^. But this is exactly the situation in which we have one massive and one 
massless mode in the theory and the low energy physics is that of a massless boson. This 
matching of the two results at large would imply that the space of retained states must 
have a non-vanishing overlap with the true low lying states of the theory for finite values of 
which is all that is needed to show that the CORE method must work. 
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FIG. 2. The energies of one-particle states (in units of 1/a) for modified SLAG derivative for 

two different values of ji. 
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